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Abstract. - In view of a recent controversy we investigated the Mott-Hubbard transition 
in D = oo with a novel cluster approach, i) We show that any truncated Bethe lattice of 
order n can be mapped exactly to a finite Hubbard-like cluster, ii) We evaluate the self-energy 
numerically for n = 0,1,2 and compare with a series of self-consistent equation-of-motion 
solutions, hi) We find the gap to open continously at the critical U c ~ 2.5t* (t = t* /v4d). 
iv) A low-energy theory for the Mott-Hubbard transition is developed and relations between 
critical exponents are presented. 



Introduction. - The Mott-Hubbard (MH) transition, as the metal-insulator transition 
in translationally invariant systems of interacting electrons is called, is little understood at 
present. The nature of the order parameter remains unresolved and a number of phenomeno- 
logical scenarios have been proposed to describe the underlying physics. In the Brinkmann- 
Rice [1] scenario the number of charge carriers drives the MH transition, another view [2] 
suggests a binding/unbinding transition of doubly-occupied and empty sites as the appropri- 
ate model. In this second picture the MH transition coincides with the convergence radius of 
a large-interaction expansion [3] . 

The study of interacting electrons in the limit of high dimensions [4, 5] has proven very 
useful, as both analytical and numerical methods simplify in the limit of infinite dimensions. 
For instance, Monte Carlo studies [6, 7] are not limited by finite-cluster effects, only by the 
imaginary-time resolution in infinite dimensions. Here we focus on the zero-temperature half- 
filled Hubbard model on the infinite-dimensional Bethe lattice in the paramagnetic state( 1 ). 
A Mott-Hubbard transition is known [7, 8] to occur as a function of interaction strength and 
has been examined by Quantum Monte Carlo [7], by the iterative perturbation theory [7-10], 
by self-consistent diagonalization studies [11], by a non-crossing approximation [12] and a 
modified-equation-of-motion approach [13]. In a recently predicted [10,11] scenario for the 
Mott-Hubbard transition the Fermi-liquid effective mass would diverge on the metallic side 
of the transition, while no precursor of the transition would be seen on the insulating side for 
any observable, e.g. that the gap would close discontinuously from a finite value to zero at 
the transition point. Here we want to examine this unusual scenario with a novel approach. 

( 1 )Note that the Mott-Hubbard transition occurs at zero temperature only in the paramagnetic sector, while 
the true ground state is antiferromagnetically ordered at half-filling. 
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Fig. 1 - Illustration of the Hubbard star, S(l), and of the star of the stars, S(2). The lines connect 
sites connected by hopping matrix elements. 

Up to date it has not been clear how to formulate a systematic exact-diagonalization 
approach in infinite dimensions, as clusters of physical interest contain an infinite number of 
sites. The smallest cluster (see fig. 1) has been shown by van Dongen et al. [14], using an 
equation-of-motion approach, to be equivalent to a 3-site Hubbard cluster. Here we show, 
for the first time, that any truncated Bcthe lattice may be mapped to an effective Hubbard- 
like cluster. The exact diagonalization of these clusters allows to calculate the local Green's 
function. We then estimate the coefficients of the Laurent expansion of the self-energy and 
show that a simple phenomenological low-energy theory for the Mott-Hubbard transition can 
be formulated in terms of a truncated Laurent expansion of the self-energy. 

Star of the stars- Any truncated Bethe lattice in infinite dimensions can be mapped 
exactly to a Hubbard-like cluster with a finite number of sites. The ground state of this 
cluster can be used to obtain the desired Green's function of the truncated lattice. By u S(n)" 
we denote a truncated Bethe lattice, which contains all sites which are separated by at most 
"n" steps from the central site. 5(0) denotes the isolated atom, 5(1) an atom connected to 
2D n.n. sites and 5(2) an atom with 2D n.n. and 2D(2D — 1) n.n.n. sites (see fig. 1 for 
an illustration). For the Hubbard- model there is an on-site repulsion U on every site and a 
hopping [4] t = t* /V4D matrix element between n.n. sites. 

We now consider the zero-temperature, retarded Green's functions G^^lo) and G^{ui) 
of the central site of S(n) and a n.n. site respectively. In the limit D — > oo all diagrams con- 
tributing to (u>) which contain the central site are suppressed as t 2 <~ 1/D — ► 0. Therefore 
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G^ ' (u>) — G " ^(u>) and the quantity of interest, G^ (oj) can be calculated recursively. 

JVn-i 

Denoting by N n the number of poles of Go"^(w) we may write Gq" ^(oj) = fj/i 1 ^ + 

j 

(i — sj), with fj > and fj = 1- ^o"" 1 ^^) enters the equation of motion for Gq 1 \lo) 

3 

always as t 2 G^ 1 Therefore S(n) is equivalent to an effective system where every one of 

the 2D n.n. sites is replaced by N n _\ non-interacting orbitals which are coupled by hopping 
matrix elements tj — ty/f] to the central site. The central site Green's function couples 
only to the respective symmetric combination of these iV„_i orbitals, as they are mutually 
non-interacting. Gq™^ (lu) is, therefore, determined by solving an effective Hubbard cluster 
where one central site (with finite U) is connected to iV n _i other sites (with on-site energies 
Sj and zero U) by hopping matrix elements = V2Dtj = t*y/ fj/2. For N — 2 we find 

G { q } (lu) = l/2[l/(w - U/2) + + U/2)] at half-filling and, therefore, 5(1) is equivalent to 
an effective 3-site Hubbard-cluster [14]. The 22 poles of the Hubbard-star Green's function, 
Gq 1 ^ (to), are reduced by symmetries to an effective N\ — 14. Using standard techniques [16], we 
have been able to obtain Gq ' (uj) at half-filling from the ground-state of the 15-site Hubbard- 
likc cluster in the parametric sector, which was obtained by a direct-diagonalization technique. 

Self-energy. - The Laurent series of the self-energy has the following functional form: 

E( W ) = a/ U + M + (1 - 1/z ) u> l _ {l _ l/z , )ujGmc{uu y (D 

where a and z* are parameters depending on U. Gi nc (u>) is an "incoherent" Green's function 
yet to be determined. Equation (1) turns out to be useful in the discussion of the low-energy 
properties of the self-energy. The Fermi-liquid state corresponds to a = 0, the MH insulating 
state to a > [8,17]. 

The first two terms of the r.h.s. of cq. (1) arise from the definition of self-energy 

V ' " 2 K > G («) (w) ' v ' 

on clusters S(n) [18] and particle-hole symmetry at half-filling which demands G^\uj) = 

N n /2 N n /2 

J2 2f 3 uo/{uj 2 - efj and, therefore, a= J2 £j/{ 2 fj) > °- 

3 3 

The third term of the r.h.s. of eq. (1) is a consequence of the fact that in infinite 
dimensions the self-energy is site diagonal. In the metallic state (a = 0) we write the Green's 
function at the Fermi level as [19] 

^ + (1 - Z)G (inc) (u;) = — . (3 ) 

Solving eq. (3 ) for ~^2(lo) as a function of Ga nc -\(u) and Z, we obtain eq. (1) with a = 
and z* = Z. G( inc )(u) is defined through eq. (1) for the case a > 0. Our numerical results 
indicate that this G( inc )(o;) has indeed all the properties of a Green's function, in particular 
that ^G(lu) < and - / oIlu^G(lu) /tt = 1. 



4 



EUROPHYSICS LETTERS 




Fig. 2 - Results from exact diagonalization of the atom (S(0), triangles), the star (S(l), open circles) 
and the star of the stars (5(2), filled circles) as a function of U in units of t* . a) The smallest pole of 
Go™'(o»), which corresponds to half the gap. Denoted by 'rec' is the prediction for the thermodynamic 

limit by the recursion method analysis, b), c) The coefficients a (a^ = o^ 3 ', a' 4 ' and 

2* (insert) of the Laurent-series of the self-energy: £(w) = a/co + U/2 + (1 — l/z*)u> + ... as given by 
eq. (1). Only for a = is z* = Z = rn/m*. The a (n) and Z (n) (n = 2, 3, 4) are the predictions of the 
equation-of-motion solutions. A 5(1), o 5(2), • 5(3). The dotted lines are guides to the eye. 



Results - We have exactly diagonalizcd the half- filled Hubbard- like Hamiltonian on S(n) 
for n = 0, 1 (for 5(2) we have 8— | and 7— J, particles corresponding to ~ 4 • 10 7 different con- 
figurations) and determined the one-particle Green's function, G^\uj), in the paramagnetic 
(see footnote (*)) superposition of the degenerate spin doublets via a continued- fraction ex- 
pansion [16]. 
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In fig. 2a) we have plotted the position of the lowest pole in Gq (uj) (as a function of U), 
corresponding to half of the optical gap, A. In order to estimate the influence of the finite 
system size on A we have subjected the Green's function for 5(2) to a continued-fraction 
analysis which has been successfully used [20] to eliminate finite-size effects in numerical 
calculations of dynamical correlation functions from small systems. The first few coefficients 
of the continued-fraction representation of the Green's function allow for reliable estimates of 
its low-energy behavior in the thermodynamic limit. We find that the gap vanishes identically 
for U < 2 in the thermodynamic limit and that it is nicely linear for U > 3, with the same 
slope as the cluster data for 5(2). These results, see fig. 2a), indicate that our estimates of 
the self-energy parameters from 5(2) are well controlled and that we may reliably estimate 
U c ~ 2M*. 

In fig. 2b) we present the results for a, as defined by eq. (1). a(U) increases between U = 2 
and U = 3 by a factor of twelve for S(2), indicative of a gap opening in the thermodynamic 
limit, 5(n — ► oo). 

In fig. 2c) we show the results for z*, as defined by eq. (1). Note that z* corresponds to 
the quasi-particle renormalization factor Z = m/m* in the thermodynamic limit whenever 
a = 0. We find for both S(l) and S(2) the minimum of z* occurs at approximatively the same 
Umin ~ 2.5t* and that z* is a continuous function of U. 

In fig. 2 we have included the results of a recently proposed systematic series of self- 
consistent equation-of-motion solutions [21], and Z^> (with n — 2,3,4). The n-th-order 
equation-of-motion solution satisfies the first n — 1 equations of motion exactly, while the n-th- 
order equation of motion is decoupled self-consistcntly [21]. n — 2 corresponds to the Hubbard- 
III solution [22]. The critical [7 C (2) = [7 C (3) = VoF - 2.45t* and [7 C (4) = ^40l* - 2.52t* found 
by the equation-of-motion solutions [21] agree very well with the U c ~ 2, 5t* estimated from 
the 5(2) cluster with the recursion method, see fig. 2a) and b), and agrees also well with 
the location of the minimum in z*(U), as obtained for both 5(1) and 5(2), see fig. 2c). The 
crowding of lines in fig. 2b) for U > 3 puts into evidence the remarkable consistency between 
the cluster results and the equation-of-motion results for a. 

In Fig. 3 we plot the results for —ImGo(u> + iS) and +ImGi nc (uj + iS) as obtained for 
the star of the stars for 5 = 0.05i* and U = 2. Gq{w) has, as expected, poles at very small 
frequencies, an indication of a gapless density of states in the thermodynamic limit. Gi nc (u)), 
on the other hand, has no poles at small frequencies, the structure in Gi nc (u>) seems to start 
at about [7/2. In the metallic state the quasi-particle damping if given by <~ oj 2 ImGi nc (0) 
and ImGi nc (0) is finite in the thermodynamic limit. In the insulating state for both Go (a;) 
and Gi nc (u) a gap A opens and ImGi nc {Q) = 0. 

Discussion - Our cluster results for A, a and z* (see fig. 2), the analysis of the cluster 
Green's function with the recursion method [20] and the series of self-consistent equation-of- 
motion solutions [21] are all consistent with i) a U c ~ 2.5t*, ii) the gap A opening continuously 
at U c and iii) the parameter z* of the self-energy representation equation (1) being a continu- 
ous function of U at the Mott-Hubbard transition. These three results do not agree with other 
studies [8]- [11], based on the iterative perturbation theory and self-consistent diagonalization 
studies, which find a) a U c — 4.2i* — 4.7t*, b) the gap A to open discontinuously at U c and c) 
the parameter z* to be discontinuous at U c . 

Effective theory - The results for the Gi nc {u>) presented in fig. 3 indicate that Gi nc {oj) is 
a smooth function, for small frequencies, in the thermodynamic limit. The first few terms of 
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a small-frequency expansion for the self-energy (compare Eq. (1)) 

i>>^ + !+H)—H)^> w 

where 7 = G, nc (0), are then sufficient to describe the asymptotic, small- frequency behavior 
near the Mott-Hubbard transition. The three parameters x = a,z*,-f entering eq. (4) are 
expected to scale like \U — U c \ Xl in the critical region \U — U c \ << U c of the Mott-Hubbard 
transition. The exact values of the critical exponents, A x , are not yet known. The equation- 
of-motion solutions [21] predict, for the Bethe lattice A^ = 1 and A 2 » = 2, the iterative 
perturbation theory and self-consistent diagonalization studies [8]- [11] A z » = 2. 

From the small- frequency expansion of the self-energy, eq. (4), one can determine the low- 
frequency behavior of the one-particle Green's function [to — (sp — fJ,) — ^(w)] -1 , where the ep 
are the one-particle eigenenergies of the Bethe lattice. For the position of the quasi-particle 
poles one finds z*e^/2[l ± </l + Aa/ (z*e^)]. For a > a Mott-Hubbard gap opens, which 

scales for A* > A a as ~ y/az* ~ \U — f7 c |( A ^+ A z*)/ 2 . Defining the scaling exponent of the gap, 
A, via A ^ |J7 — U C \^Q(U — U c )2 we then arrive at the prediction 

2A A = A, + \ z , . (5) 

We then find for the equation-of-motion solutions Aa =3/2 while the iterative perturbation 
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theory and self-consistent diagonalization studies [8]- [11] predict Aa = 0. 

Conclusions - We have shown that the properties of the Hubbard star of the stars can 
be obtained by exact diagonalization of certain 15-site Hubbard clusters. We have presented 
estimates for the optical gap A, and the quasi-particle spectral weight z*. We find our estimate 
of U c ~ 2.5t* to be in good agreement with the predictions of a recently proposed series of 
self-consistent equation-of-motion solutions [21]. 

We used the results obtained from the cluster calculations to estimate qualitatively the 
low-energy behavior of the self-energy which in turn led us to propose a simple, phenomeno- 
logical theory for the small-frequency behavior of the Mott-Hubbard transition, in terms of a 
truncated Laurent expansion of the self-energy. 
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